Asymmetry gap in the electronic band structure of bilayer graphene 
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A tight binding model is used to calculate the band structure of bilayer graphene in the presence 
of a potential difference between the layers that opens a gap A between the conduction and valence 
bands. In particular, a self consistent Hartree approximation is used to describe imperfect screening 
of an external gate, employed primarily to control the density n of electrons on the bilayer, resulting 
in a potential difference between the layers and a density dependent gap A(n). We discuss the 
influence of a finite asymmetry gap A(0) at zero excess density, caused by the screening of an 
additional transverse electric field, on observations of the quantum Hall effect. 

PACS numbers: 73.63.-b, 71.70.Di, 73.43.Cd, 81.05.Uw 



Recently there has been a huge amount of theoretical 
interest in the low energy electronic properties of ultra- 
thin graphite films, including graphene monolayers [J 
and bilayers [!, 0] . This activity followed the successful 
fabrication of ultrathin graphite films [B[ and subsequent 
measurements of an unusual sequencing of quantum Hall 
effect plateaus [f| @] • While the Hamiltonian for low en- 
ergy electrons in a monolayer describes Dirac-like chiral 
quasiparticles with Berry phase 7r [H, 0| , the low energy 
Hamiltonian for a bilayer describes chiral quasiparticles 
with a parabolic dispersion and Berry phase 2ir ,3j. 

Theoretical Q and experimental Q studies of bilayer 
graphene have shown that asymmetry A = e 2 — e± be- 
tween on-site energies €\, €i in the layers leads to a gap 
between the otherwise degenerate conduction and valence 
bands as shown schematically in Fig. [IJa). Moreover, 
in contrast to graphene monolayers, there is an exper- 
imental possibility of controlling the magnitude of the 
gap A in the spectrum of bilayer graphene through the 
use of an external gate that is employed primarily to 
control the density of electrons n on the bilayer system 
[E B S HI- Such a dependence A(n) owes its existence 
to the Coulomb interaction between electrons, but the 
density distribution itself is dependent on the value of 
A via the band structure. Here we use a self-consistent 
Hartree approximation to determine the electronic dis- 
tribution on a bilayer in the presence of an external gate. 
The numerically calculated A(n) is shown in Fig. [TJb) 
(solid line) for typical experimental parameters in the 
case that the gap A = at zero excess density [lO| |. It 
shows that the addition of density n ~ 10 12 cm~ 2 yields a 
gap A ~ lOmeV. The numerical curve is compared with 
the following analytic approximation (dashed line) valid 
for A < 71: 
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where Cb is the capacitance of a bilayer of area 
v is the in-plane velocity of the bilayer and 71 is 
the inter-plane coupling. The dimensionless parameter 
A = e 2 L 2 7i/ (2irh 2 v 2 Cb) describes the effectiveness of the 
screening of the bilayer. The limit A — > describes poor 
screening when the density on each layer is equal to n/2 
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FIG. 1: (a) schematic of the electronic bands near the K point 
in the presence of finite layer asymmetry A (for illustrative 
purposes a very large asymmetry A = 71 is used), (b) nu- 
merically calculated dependence A(n) (solid line) compared 
with analytic expression (dashed line) using Eq. |T} . For these 
values of density A w A fiol ]. 



whereas for A — > 00 there is excellent screening, the den- 
sity lies solely on the layer closest to the extern al g ate 
and A — 0. For typical experimental parameters 11| we 
estimate that A ~ 1. 

We consider a graphene bilayer, with interlayer sepa- 
ration Co, located a distance d from a gate. An external 
gate voltage V g = end/e' r eo induces a total excess density 
n = ni + ri2 on the bilayer where ri\ (n^) is the excess 
density on the layer closest to (furthest from) the gate. 
Here £0 is the permittivity of free space, e' r is the dielec- 
tric constant and e is the electronic charge. Imperfect 
screening of the effective charge density en from the gate 
leads to an excess density rii on the layer furthest from 
the gate, with a corresponding change in potential energy 



AU2 = e ri2Co/e r eo that determines the layer asymmetry 



A (n) = e 2 - ei = A + e 2 n 2 L 2 /C b 



(2) 



where Cb = e r EoL 2 /cq and e r is the bilayer dielectric 
constant. We introduce the bare asymmetry parame- 
ter Aq due to an additional transverse electric field pro- 
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ducing finite asymmetry A(0) at zero density. The fol- 
lowing analysis includes (a) the tight binding model of 
bilayer graphene including the asymmetry A, (b) a self- 
consistent determination of the gap numerically and, for 
£f,7i 3> A, analytically, and (c) a description of the in- 
fluence of Ao on the cyclotron mass and the sequencing 
of quantum Hall effect plateaus at low density. 

(a) The bilayer is modelled as two coupled hexagonal 
lattices with inequivalent sites Al, Bl and A2, B2 in the 
first and second graphene sheets, respectively, arranged 
according to Bernal (A2-B1) stacking. This lattice has 
a degeneracy point at each of two inequivalent corners, 
K and K, of the hexagonal Brillouin zone which deter- 
mine the centers of two valleys of a gapless spectrum. 
At the degeneracy point, electron states on inequivalent 
(Al/Bl or A2/B2) sublattices in a single layer are de- 
coupled, whereas interlayer coupling JA2B1 = 7i forms 
dimers from pairs of A2-B1 orbitals in a bilayer, thus 
leading to the formation of high energy bands [U, ll2T |. 
The tight-binding Hamiltonian operates in the space of 
wave functions VP = (ipAi,ipB2, ipA2,ipBi) in the valley K 
and of — (^B2,'<pAi>'^Bi,ipA2) in the valley K @: 
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where n = p x + ip y and £ = +1(— 1) labels valley K (K). 
The Hamiltonian takes into account asymmetry A = e 2 — 
ei between on-site energies in the two layers, e 2 = 5A, 
€1 = — i A. Within each plane, nearest neighbor coupling 
as parameterised by 70 leads to in-plane velocity v = 
(V3/2) a7o/fi. where a is the lattice constant. 

At zero magnetic field, the Hamiltonian Ti describes 
four valley-degenerate bands, e± (p), a = 1,2, with 



e(Q)2 = y + t + vV + (_ir V t + y2p2 (7? + A2) ' 

where p is the magnitude of the momentum near the K 
point, Fig. [Ha). The energies of the bands exactly at 
the K point are |e^ 2) (0)| = ^7? + A 7 4 and je^O)] = 
|A|/2. Thus ej_ describes higher-energy (A2-B1 dimer) 

bands |e^'| > 71 whereas e± are low ener gy b ands split 
by the layer asymmetry A at the K point [10| . 

The electronic densities n\ and 77,2 on the in- 
dividual layers are given by an integral with re- 
spect to momentum over the circular Fermi sur- 
face n H2) = (2/irh 2 ) Jpdp(\iP A1{2) {p)\ 2 + |i/>bi ( 2)(p)| 2 ) 
where we have included a factor of four to take into ac- 
count spin and valley degeneracy. On determining the 
wavefunction amplitudes on the four atomic sites we find 
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where the minus (plus) sign is for the first (second) layer. 
The limits of integration are allowed values of momentum 
that, depending on the band in question and the value 
of the Fermi energy tp, are p = or p± where v 2 pj n — 
4 + iA 2 ±[4(7 1 2 + A 2 )-l7 1 2 A 2 ] 1 /2. 

(b) The self-consistent calculation begins with zero ex- 
ternal gate voltage, assuming the Fermi energy to be 
located at e = 0. We take into account the contribu- 
tion of the two valence bands to the initial densities on 
each layer, nf^ = —n 2 , and evaluate the gap A(0) us- 
ing Eqs. (|2I4|) . Then, we proceed to the case where an 
external gate voltage produces excess densities m, n 2 , 
n = m + n 2 7^ 0, and a gap A(n). Figure a) shows 
A(n) for different values of the bare asymmetry Ao 
and Fig. 12b) shows the variation of the individual layer 
densities n± and n% for finite Ao. The main effect of finite 
A (dashed and dotted lines in Fig. [2]^a)) is to shift the 
plot A(n) with respect to the Ao = case (solid line), 
but, as a result of screening, the asymmetry at zero den- 
sity A(0) is smaller than the bare asymmetry Ao. For ex- 
ample, the dotted curve for Ao = O.271 = 78meV shows 
A(0) » 21meV. Figure [He) shows the variation of A(0) 
with bare asymmetry Ao (solid line). 

We find an analytic expression for A(n) in the limit 
71, ep |A| by evaluating the density on each layer 
Eq. ((4]) using an expansion in A/71. For simplicity, we 
consider |A|/2 < |ejr| < 71 so that either the band just 
above or just below zero energy is partially occupied. On 
integrating Eq. Q from zero momentum to pp — p+, we 
find the densities in the partially occupied bands. We also 
take into account the redistribution of density within the 
valence bands for ep ^ as compared to the ep = case, 
as given by ±71 A/ (4irH 2 v 2 ) In (471/I A|), so that 
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where eo = (71/2) [yT + <iv 2 p' F /"/ 2 — l] and n = ni+n 2 ~ 
sgn(e F )p 2 F / {irh 2 ). If |A| < j|e F |, 7 i} so that v 2 p 2 F w 
ep + 7i|c_f I j then eo is approximately independent of A, 
and we may use Eq. @ to find the gap, 
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1 + A [(eo/7i) + (eo/7i) 2 - I In (eo/7i)] 
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This is valid for intermediate densities |A|,|A | -C \ep\ < 

71 where ep ~ ±eo [at very low density it incorrectly pre- 
dicts A(0) = for A 7^ 0]. In the regime AitK 2 v 2 \n\ < 

7 2 and A = it gives Eq. (TTJ). 

There is uncertainty in the value of the capacitance 
Cb within this model because we assume that the ex- 
cess charge is uniformly distributed within infinitesimally 
thin 2d layers. Also, we used s r = 1 for the dielec- 
tric constant of the bilayer which agrees with the pre- 
diction for other low dimensional structures [l3|, but is 
smaller than the value for bulk graphite e r « 2.4 
Furthermore, we assume that the model parameters, in- 
cluding layer separation and interlayer coupling 71, do 
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FIG. 2: Numerical evaluation of (a) the bilayer asymmetry 
A(n) for different values of the bare asymmetry Ao = (solid 
line), Ao = O.I71 = 39meV (dashed line) and Ao = O.271 = 
78meV (dotted line), using typical parameter values [llj], (b) 
layer densities ni (solid) and n\ (dashed) as a function of n 
for A = O.271 = 78meV, (c) A(0) as a function of A for 
£r = 1 (solid line) and e r = 2 (dashed), and (d) the cyclotron 
mass m c in units of the bare electronic mass m e , for different 
values of Ao as in (a). 



not vary as a function of applied potential whereas a re- 
cent experiment j8J found that 71 changed by about 3% 
as A changed by lOOmcV. The effectiveness of screening 
A = e 2 L 2 7i/(27rfi 2 z; 2 C{,) depends on the model param- 
eters as illustrated in Fig. [2jc) which shows A(0) as a 
function of Ao for e r — 1 (solid line) and e r — 2 (dashed): 
the gap A(0) increases with the dielectric constant. 

(c) We use the self consistent analysis to evaluate 
the semiclassical expression for cyclotron mass m c that 
may be measured in an experimental observation at fi- 
nite magnetic field of, say, Shubnikov de Haas oscillations 
@, fl 0- For a circular Fermi surface, the semiclassical 
expression for cyclotron mass m c — p/(de/dp) gives 
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In the limit A = for the low energy b ands a = 1 this 
gives m c = ( ll /2v 2 )y/l + Ai:h 2 v 2 \n\hl 0. Figure E^d) 
shows the cyclotron mass m c as a function of total density 
n. For Ao = (solid line) the mass is symmetric and 
finite m c = 71 /2v 2 at zero density. The other curves, for 
A > 0, show that the cyclotron mass for positive and 
negative densities is asymmetric with divergent behavior 
at low density resulting from the "mexican hat" structure 
of the low energy bands [HJ ■ 

Another manifestation of finite Ao is in the sequencing 
of quantum Hall effect (QHE) plateaus at low density. 



The Landau level (LL) spectrum of bilayer graphene has 
been discussed in Ref. [3j: here we briefly mention the 
features relevant for finite Ao. The low energy states 
of electrons in bands e£ are conveniently described by 
an effective two component Hamiltonian that operates in 
the space of wave functions W = (tpAi, ipB2) in the valley 
K and of \& = (V>B2, "4>Ai) in the valley K [3[: 
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where n = p x + ip y and B is the magnetic field. For com- 
pleteness we include the final two terms, with dimension- 
less parameters \a\, |/?| <C 1, describing weak asymmetry 
between the conduction and valence bands. Neglected in 
Ref. 3, they arise due to additional couplings such as 
next-to-nearest neighbor in-plane coupling between ^41- 
Al and B2-B2 orbitals or intcrlayer coupling between 
A1-A2 and B1-B2 orbitals. The Landau level (LL) spec- 
trum may be found using the Landau gauge A = (0, Bx), 
in which operators tp and it coincide with raising and 
lowering operators [HI in the basis of Landau functions 



N 



•n{x), such that ir^ 
= -i(H/\ B )V2Nct>N 



N = i{h/\ B )y/2(N + l)<l> N+1 , 
and tt 4>o — 0, where Xb = 
yJh/{eB). The spectrum has almost equidistant energy 
levels for N > 2 with spacing hu> c , lo c = eB/m, which 
are weakly split in valleys K and K (£ = ±1) 0], 
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where and are assigned to electron and hole states, 
respectively, and we assume |A|/7i <C 1, \A\/fku c < 1. 

The LL spectrum in each valley also contains two levels 
identified using the fact that ir 2 <fii = tt 2 (/>q = 0, 



= -^A 
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(^0,0); 
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According to different definitions of two-component <& in 
two valleys, N = 0,1 LL states in the valley K are formed 
by orbitals predominantly on the A sites from the bottom 
layer, whereas the corresponding states in the valley K 
are located on B sites from the top layer, as reflected by 
the splitting A between the lowest LL in the two valleys. 
The spectrum is shown on the right side of Fig. [3] assum- 
ing |A|/7i -c 1, \A\/Huj c < 1. In a symmetric bilayer 
(A = 0) levels e and e\ are degenerate and have the 
same energy in valleys K and K, thus forming an 8-fold 
degenerate LL at e = (here, spin is taken into account) 
as shown on the left side of Fig. [3] 

The LL spectrum is reflected in the QHE Hall conduc- 
tivity dependence on carrier density, a xy (n), although, 
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FIG. 3: Landau levels for zero (left) and finite asymmetry 
A (right). Brackets (iV, £) indicate LL number N and valley 
index £ = ±1. In the center the predicted Hall conductivity 
a xy as a function of carrier density for bare asymmetry Ao = 
(dashed line) is compared to that for Ao / (solid line). 



in general, finite temperature and LL broadening tend to 
mask small LL splitting. When the asymmetry at zero 
density A(0), Fig.[2^c), is large enough that the percolat- 
ing states [15| from the lowest levels in the two valleys can 
be resolved (the right side of Fig. [3]), then these levels, 
embedded into the ladder of 4- fold degenerate LL's with 
n > 2, Eq. result in the form of o xy (n) as sketched 
in the center of Fig. [3] (solid line) which exhibits plateaus 
at all integer values of 4e 2 /h including a plateau at zero 
density. The temperature dependence of the zero density 
plateau will differ from that of the other plateaus because 
of the different activation energies, related to A(0) and 
hu) c , respectively. When it is not possible to resolve the 
asymmetry splitting of the group of 8 states at |e| =0 (for 
A(0) « A « 0) then a xy (n) (dashed line in Fig. [3]) has 
plateaus at integer values of 4e 2 /h and a double 8e 2 /h 



step across n = d, [ 

The use of doping [8j recently achieved gaps of up to 
about 200meV. Alternatively, a potential difference of 
30V applied to a back gate placed 300nm from the bilayer 
(combined with — 3V applied to a top gate 30nm away) 
translates into 30mV across the bilayer, Ao = 30mcV, 
and A(0) ss 7meV. This compares with Huj c S3 22meV for 
field B = 10T and level broadening due to scattering 
of approximately 3meV in recent experiments @ . 

This paper has considered the influence of a single gate 
placed above an ultra thin graphite film as employed in 
recent experiments [!, H, 0,0 to control the density of 
electrons n on the film. For the typical experimental 
range of densities shown in Fig. \V[b) the dependence of 
the asymmetry gap A(n) is roughly linear with n with 
the addition of density n ~ 10 12 cm~ 2 yielding a gap 
A ~ lOmeV. Finite bilayer asymmetry gap A(0) at zero 
excess density, due to an additional transverse electric 
field, could influence observations of the integer quan- 
tum Hall effect by introducing a plateau at zero density 
in the Hall conductivity. This plateau is accompanied 
by a dip in the diagonal conductivity, the temperature 
dependence of which can be used to extract the value of 
the asymmetry gap A(0). This is distinguishable from 
the behavior of the higher plateaus with activation ener- 
gies determined by the cyclotron frequency. 
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